Ferromagnetic superconductivity driven by changing Fermi surface topology 
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We introduce a simple but powerful zero temperature Stoner model to explain the unusual phase 
diagram of the ferromagnetic superconductor, UGe2. Triplet superconductivity is driven in the ferro- 
magnetic phase by tuning the majority spin Fermi level through one of two peaks in the paramagnetic 
density of states (DOS). Each peak is associated with a metamagnetic jump in magnetisation. The 
twin peak DOS may be derived from a tight-binding, quasi-one-dimensional bandstructure, inspired 
by previous bandstructure calculations. 



I. INTRODUCTION 

The predominance of spin singlet superconductors over 
their triplet counterparts has, in part, lead to the be- 
lief that superconductivity and magnetism are mutually 
exclusive — the upper critical field of a singlet supercon- 
ductor is bounded by the Pauli paramagnetic limit. One 
might expect that a ferromagnet would be the natu- 
ral stage for spin triplet superconductivity, where we 
could overcome Pauli limiting effects. But, until recently, 
there were no examples of 'ferromagnetic superconduc- 
tivity' (FMSC) — the coexistence of itinerant ferromag- 
netisrn (FM) and superconductivity (SC) in a single bulk 
phasea. This situation has changed with, the observation 
of FMSC in UGe 2 B, URhGeB and ZrZna^. The behaviour 
of these materials is an example of a more general phe- 
nomenon; the observation of a novel state on the bor- 
der of magnetism at low temperatures. By supressing 
magnetic order bj£ some control parameter, be it elec- 
tron/hole densitya or external pressures physicists are 
now able to access regimes where magnetic fluctuations 
become quantum critical. 

We choose here to concentrate on the case of UGe2, 
because whilst SC is only measurable in the ferromag- 
netic state (in common with URhGe and ZrZn2), UGe2 
seems to posess particularly low electronic dimensional- 
ity, uniaxial magnetisation, and revealing features in the 
temperature-pressure phase diagram which we now re- 
view. 

In Fig. [I] we show the temperature-pressure phase di- 
agram for UGe2, with the Curie temperature Tq (su- 
pressed to zero at pressure p c ) ancL»Aoerconducting tran- 
sition temperature T$c indicatecootl. Another feature, 
T x is also shown— This T x shows up in measurements of 
lattice expansionEl, as a jump in the low temperature T 2 
component of the resistivityn, as a small enhancement in 
specific heaifl, as a kink in resistivitya and as a change 
in the character of the Fermi surface as measured in de 
Haas van Alphen experiments^!]. Most importantly for 
thisrivork, T x also appears as a slight jump in magnetisa- 
tionfcl which is sharpened at lower temperatures such that 
the low temperature moment has a step at pressure p x 
in addition to the step at the quantum phase transition 
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FIG. 1: The temperature-pressure phase diagram a uU and low 
temperature easy axis magnetic moment, M (after Pfleiderer 
and HuxleyB) of UGe2. Tsc is the superconducting transition 
temperature, scaled by a factor of 10 for clarity. Tc denotes 
the Curie temperature. T x is a feature in the ferromagnetic 
state seen in various measurements, described in the text. M 
shows two transitions; one at p x , corresponding to where the 
T x line meets the pressure axis, and the other at the critical 
Curie pressure, p c . 



pressure, p c (see Fig. Q). Furthermore, we note the close 
proximity of the T x (p) line to the peak in Tsc- 

Most theories which describe ferromagnets close to a 
quantum phase transition have predicted that the su- 
perconducting transition temperature, Tsc should be at 
least as high in the paramagnetic state as it is in the fer- 
romagnetic state. These theories have considered a elec- 
tronically three^dimensionaLfprromagnet, either magnet- 
ically isotropicE3 or uniaxials. Where theoretical models 
have predicted an enhancement of Tsc in the ferromag- 
netic regime, their basis seems unjustified in the case 
of UGe2. Kirkpatrick and coworkeralJ have predicted 
an enhancement of the superconducting Tsc due to the 
coupling of magnons to the longitudinal magnetic sus- 
ceptibility. However, the ferromagnetic state of UGe2 
is so magnetically anisotropic that the presence of trans- 
verse magnons seems an unlikely primary explanation for 
the enhancement of Tsc — at 4.2K and an external mag- 
netic field of 4T, the magnetisation along the easy axis 
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is about 20 to 30 times that along either of the other 
crystallographic axesli3. Other authors have drawn their 
inspiration frpm-h&ndstructure calculations. Bandstruc- 
ture analysesEj'LrEj of UGe2 seem to indicate that low- 
ering temperature sparks the evolution of a quasi-two- 
dimensional majority carrier Fermi surface sheet below 
Tc- Furthermore there is the possibility that large sec- 
tions of the quasi-two-dimensional Fermi surface may be 
parallel, making it almost one-dimensional. Until now, 
this low-dimensional magnetism has pushed authors in 
the direction of postulating the existence of a chaEgg- 
or spin-density-wave state (CDW or SDW) below T x an&, 
sometimes by analogv— .with the a phase of uranium. 
Watanabe and Miyakal3 have postulated that the inter- 
play of CD or SD fluctuations at high wavevector will 
couple to the magnetisation, M in such a way as to en- 
hance it at some critical value, MjUl. However, spin 
density fluctuations have yet to be observed in neutron 
experiments on UGe2- 

We turn to the low dimensional bandstructure for a dif- 
ferent effect. The key idea will be that in a ferromagnet, 
somewhat uniquely, the magnetisation acts as a tuning 
parameter which can subtly change the topology of the 
anisotropic Fermi surfaces of different spin species. By 
contrast, in a rigid band picture of a paramagnetic metal, 
the Fermi surface is fixed. The added topological possi- 
bilities for a ferromagnet should be viewed as a useful 
tool — and as a reason for observing the enhancement of 
features, such as Tsc, within the ferromagnetic phase. 
This paper is planned as follows: firstly, we show that an 
electronic density of states (DOS) which has two peaks 
can reproduce the two steps in the observed low temper- 
ature magnetisation. We then show that the necessary 
form of DOS arises naturally from a low (quasi-one-) di- 
mensional bandstructure and that the magnetisation re- 
sulting from this bandstructure has a jump in the ferro- 
magnetic state which is coincident with the maximum in 
a superconducting instability, mediated by spin fluctua- 
tions. 



II. MODEL 



ever the reduced electronic dimensionality of UGe2 will 
probably heighten the importance of spin fluctuations at 
finite temperatures. 

We therefore circumvent finite temperature concerns 
by employing a zero temperature Stoner theory, with the 
first aim being to reproduce the step in M(p) at p x . We 
consider the action of pressure to be akin to that of vary- 
ing the exchange energy, / in a Stoner model of the one- 
electron energy of separated majority (say, |) and minor- 
ity (say, |) spin sheets: E ka = e k ± IM (- f, + I)- In 
this description, the occupation of each spin sheet a is 
n a = p(e)de, where M = |(n-f — n\). Here €b is the 
bottom of the band, and spin a occupies energy states up 
to fj, a . The DOS is given by p(e). We consider the total 
number of spins, rij + = N to be fixed. The chemical 
potential, p a of each spin sheet is therefore completely 
determined by the particular magnetisation and chosen 
electron number. The total energy density of the electron 
system is 

F[M] = / ep{e)de + / ep{e)de 
N 2 

+/( — - M ) - gp B HM, (1) 

where we have included a term for the presence of an 
external magnetic field, H. 

Most phenomenological expansions of this energy den- 
sity have included terms even in M, up to order M 6 . 
This can give one first order transition in M(I). How- 
ever, we are looking for an additional transition, corre- 
sponding to p x , and believed to be first ordera, although 
there is some controversy over thisE3. To have the pos- 
sibility of an additional first order transition, we need 
the next even term in the Landau free energy expansion, 
thus obtaining an M 8 theoryEl We now show that a 
DOS with two peaks gcncrically brings about the M 8 
term in F[M] by allowing a scenario where the global 
magnetisation can change rapidly, twice. We begin by 
taking a one-band DOS comprised of two Lorenztians, 
normalised such that the maximum number of electrons 
of each spin in the band is 1: 



The simplest Stoner theory of magnetism is long- 
known to be inadequate in describing the tempera- 
ture dependence of the magnetisation even of magnet- 
ically isotropic, electronically 3-dimensional ferromag- 
nets, where spin fluctuations have to be taken into ac- 
count. There, a fluctuation-averaged equation— of state 
method, in the spirit of Lonzarich and TailleferEJ or Ya- 
madaO would be a improved model of magnetism at fi- 
nite temperature, where we expect the temperature de- 
pendence of M to arise from the fluctuation response, 
rather than just the Fermi functions included in Stoner 
theory. The pronounced magnetic anisotropy of UGc2 
might ordinarily simplify matters, as the consequent ab- 
sence of transverse spin modes will make a Stoner ap- 
proach more valid, especially at low temperatures. How- 
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where we can vary a to adjust the sharpness of the two 
humps and 6 to adjust their position. The humps are 
symmetric with respect to the zero energy and are cen- 
tred on ±6. We minimise F[M], given by Eq. [IJwith re- 
spect to magnetisation, M to obtain the variation of M 
with respect to the Stoner exchange, /. In Figure || we 
show plots of M(I) for three sets of parameters a and 6, 
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Fermi level is off-centre with respect to the two DOS 
peaks. The main result is that, by making the density of 
states peaks sufficiently sharp and close together, we can 
obtain two magnetic transitions, one from the paramag- 
netic state and another within the ferromagnetic state. If 
the DOS peaks are too far apart, a saturated magnetic 
state is favoured; if the peaks are too smooth, the tran- 
sitions in M(7) are weakened and if the starting param- 
agnetic level is too close to the half-way point between 
the two peaks (N too close to 1) then only one transi- 
tion is observed. Therefore, to observe two transitions, 
the paramagnetic filling level should be between the DOS 
peaks, but off-centre with respect to them. 



A. Bandstructure phenomenology 




FIG. 2: Using the twin-peak density of states given by Eqs. |2j 
and [], it is possible to obtain two transitions in the magneti- 
sation of a Stoner ferromagnet, where the exchange parameter 
/ is the varied quantity. We plot both the normalised den- 
sity of states and the resultant M(I) curve. If the two peaks 
in the DOS are too far apart, a saturated magnetic state is 
favoured. If the two peaks are too smooth, the transitions in 
M(I) are weakened. Fhe bandstructure parameters used in 
each case are indicated, together with the range of number of 
spins, N, used in the calculations of M(I). 

in each case setting the bottom and top, eb and ex of the 
band to be —2 and +2 respectively. For each set of pa- 
rameters, we show the effect of varying the paramagnetic 
Fermi level, that is to say, the number of spins, N in the 
system, in the neighbourhood of half-filling. This places 
the paramagnetic Fermi level, hef between the two DOS 
peaks. On increasing the exchange term, /, the spin sea 
is polarised and the majority Fermi level moves up from 
that of the paramagnetic state and the minority level 
moves down. The onset of ferromagnetism is governed 
by the usual Stoner criterion, Ip(^ef) = 1 in the case 
of a second order onset, or by the finite-magnetisation 
Stoner criteriorO for a first order onset: 



The finite magnetisation criterion holds for all M once 
the system is polarised, irrespective of the order of the 
onset. Each Fermi level feels the effect of the density of 
states peaks at different /, as long as the paramagnetic 



Having shown that a double-humped DOS is perhaps 
key to understanding the magnetic properties of UGe2 
at low temperatures, we now set about working towards 
a tight-binding picture of the bandstructure of this com- 
pound which can reproduce both its magnetic and su- 
perconducting properties. There have been two major 
efforts towards bandstructure calculations of UGe2 . The 
first, originally on a crystal structure with an incorrect 
space group was performed by Yamagami et al.EH and 
has since been revisedlMi This produced spin-separated 
Fermi sheets in the ferromagnetic state, centred on the 
r point, with the majority sheet strongly nested. The 
second approach, already mentioned- has been an LDA 
+U method, by Shick and Pickets, using U as a fit- 
ting parameter for the zero temperature, ambient pres- 
sure magnetisation. Again, strongly nested Fermi sheets 
were obtained, but they were of mixed spin type and 
the sheet considered by those authors as being most im- 
portant for superconductivity is centred on the crystallo- 
graphic M point in their model. Furthermore, there was 
initial disagreement on the direction of the nesting vector 
Q relative to the easy axis of magnetisation a — both 
are 'in-plane', but Yamagami considered Q i a whilst 
Shick found Q || a. In summary, we take the following 
minimal, but key ingredients in our model: 

• quasi-one-dimensional bandstructure (crucial to 
what follows) 

• strong uniaxial magnetic anisotropy and 

• spin-snlit Fermi sheets (even in the LDA +U 
workli3, there is a large exchange splitting, of the 
order of 1 eV). 

In our model, we utilise the interaction potential for 
spin fluctuation mediated pairing in_the ferromagnetic 
state, as derived by Fay and Appelij. We also follow 
their sign convention, namely that an attractive potential 
between like spins is positive. The interaction potential 
is heavily dependent on the Lindhard response, Xffff(q) 
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for the bandstructure under consideration: 



i-/ 2 x^(q)xL°L.(q) 



(5) 



where / is again the repulsive Hubbard-type contact in- 
teraction acting between opposite spins and Xct^Ci) 18 
given by 

2_, e (k + q)-e(k) ' [) 

fa°\k.) being the Fermi occupation function for spin a 
which is at chemical potential fi a (fi a is defined previ- 
ously). According to our chosen interaction mechanism, 

Eq. [|, a large x^J (q) might naively be expected to lead 
to an enhanced T$c- However, the subtlety here is that 
the interaction potential in the ferromagnetic state mixes 
longitudinal susceptibilities of majority and minority spin 
sheets so the effect is not so clear-cut. We recall that in 
the free electron model, only in electronic dimensions less 
than two, there can be a peak in Xo°i (q) at non-zero q@. 
In that case, we might expect that when one sheet of the 
Fermi surface is at optimal nesting, then perhaps V aa (q) 
will be highest. The problem with this argument is that 
a V^o-(q) dominated by high-Q modes would normally 
lead to very weak triplet pairing and in tight-binding ap- 
proaches, finite q peaks in x<to- (q) are possible in two and 
three dimensions. 

The approach we take is to look for density of states- 
driven superconductivity, where the large density of 
states giving rise to the magnetisation step at M x is 
also able to enhance superconductivity in the ferromag- 
netic state. In tight-binding theory in two dimensions 
and lower, the peak in the density of states comes from 
the presence of a van Hove singularity. The simplest, fa- 
miliar tight-binding bandstructure for the cuprates is of 
the form 

e(k) = —a x cos k x — (5 cos k x cos k y — a y cos k y (7) 

with a x — a y = 1 and p less than 1. This corresponds 
to equal nearest neighbour hopping in each of the two di- 
mensions, and includes the next-nearest neighbour term, 
cos k x cos k y . This structure, as shown in Figure|3|(a), has 
one van Hove singularity and is therefore of little use in 
studying the magnetic properties of UGe2 — as we have 
already found, we require two density of states peaks. 
However, as we will now show, if we now reduce the am- 
plitude of the cos k y term, we cross over into quasi-one- 
dimensions and begin to lift the degeneracy of the van 
Hove contours, giving two of the them, with adjustable 
separation. 

The most common quasi-lD bandstructure in the lit- 
erature is probably 







(a) 



(b) 



(c) 



e(k) 



-a x cos k x 



a y cos k y 



(8) 



FIG. 3: Two dimensional and quasi-one-dimensional tight- 
binding bandstructures: contour plots (top row of figures) 
with the van Hove contours shown in red and associated densi- 
ties of states (bottom row). In (a) the two-dimensional struc- 
ture given by Eq. |^with a x = a y = 1 and (3 = —0.3 is shown. 
This is a structure most often studied in connection with the 
cuprates. In (b) we show the quasi-one-dimensional structure 
of Eq. ^| with a x = l,a y = 0.3. In (c), higher harmonics 
in the principal direction are included (Eq. H) with a x = 1, 
P = 0.7, S = 0.03 and 7 = -0.03. We see that in (a) there is 
one (degenerate) van Hove contour, but on going into quasi- 
one-dimensions, we lift this degeneracy and find two van Hove 
contours and hence two maxima in the density of states. Fur- 
thermore, the choice of bandstructure in (c) has an optimal 
nesting scenario (shown in blue) which is separate from both 
of the van Hove contours, due to the elimination of the cos k y 
term in the bandstructure. 



which is used to describe some organic superconductors. 
Here, a x is 1 and a y is less than 1 and the model corre- 
sponds to nearest neighbour hopping in each of the two 
dimensions. This bandstructure has two van Hove con- 
tours (see Figure |^(b) and hence two density of states 
peaks. Thus, on tuning our magnetisation, we can ob- 
serve in both spin types the transition through a van 
Hove singularity in the density of states at the Fermi 
level, provided, as stated before, that our paramagnetic 
Fermi level is between the two peaks in the density of 
states. The problem arises that when the peaks are close 
enough together to cause an effect in the magnetisation 
curves, the variation of the DOS in between them is not 
very rapid. This leads to a weakening of the transitions 
in M(I) as found with a slowly varying Lorentzian DOS 
in the last section. 

However, if we assume a quasi-one-dimensional disper- 
sion of the form 



e(k) = — a x cos k x — (3 cos k x cos k y — 7 cos 2k x — 6 cos 3k x 

with a x = 1 and /3,j, and S all less than unit magnitude, 
we will be able to explore the possibility of two, first order 
jumps in M{I). 
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The higher harmonics in the principal direction k x and 
the next-nearest term, cos k x cos k y which appear in Eq. 
^ are not unreasonable in a system which we know to 
be strongly one-dimensional in character. The nearest 
neighbour term in k y has been reduced to zero. As 
shown in Figure ||(c) , this dispersion relation yields two 
van Hove contours (and hence maxima in electronic den- 
sity of states), and also goes through a perfect nesting 
scenario .Eil We note that our choice of bandstructure is 
based on an extrapolation from the one point (2Z~ r 
p = 0) of the phase diagram calculated in RefsliStLZlE 
with the idea being that there should be strong nest- 
ing present at full magnetisation in our model, to match 
the bandstructure calculations. Indeed, the high energy 
contours of this dispersion are fairly well nested (Figure 
||(c)). To achieve a magnetically saturated state which 
does not completely fill the band, we require less than 
half-filling of the band in the paramagnetic state, a re- 
quirement which is in line with the demands of the M(I) 
profile from Section ||. Our choice of bandstructure will 
ultimately help to link p x , the maximum in T$c and the 
mass enhancement observed in de Haas van Alphen mea- 
surements. First we concentrate on the magnetisation as 
a function of Stoner exchange, I. 



III. RESULTS 
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A. Two jumps in magnetisation from a 
tight-binding bandstructure 

Taking the quasi-one-dimensional bandstructure given 
by Eq. ^ and varying parameters (3, 7, and 5, we obtain 
the variation of zero temperature magnetisation with ex- 
change-interaction / as we did for Lorentzian densities of 
states.cZI In Figure f|, we present M(7) for different total 
numbers of electrons, N below half-filling in the case of 
four different sets of tight-binding parameters. In each 
case, the paramagnetic Fermi level sits in between the 
two peaks in the DOS. Initial polarisation at I c (labelled 
by analogy with p c in experiment) is first order — this is 
due to the minority sheet approaching the lower van Hove 
peak as polarisation commences. 

There are several features to be noted from the plots 
in Fig. m Firstly, as in Fig. |^, if N is below half- 
filling, a second magnetic transition can occur when the 
majority spin sheet feels the effect of the upper peak 
in the DOS. However, when N is too close to half fill- 
ing, or when the dip between the two density of states 
peaks is too sharp, we get only one first order transi- 
tion, which can even be straight to saturation magneti- 
sation. Once again, in order to observe the two tran- 
sitions seen in UGe2, one should have a paramagnetic 
Fermi level which is off-centre with respect to two DOS 
peaks. Secondly, the upper transition (which we hence- 
forth label I x by analogy with p x from experiment), can 
either be to finite magnetisation (and either first order or 
beyond the first order critical point) or can be first order 
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FIG. 4: A set of phenomenological quasi-one-dimensional dis- 
persions for UGe2, where both the possibility of tuning the 
magnetisation through a van Hove singularity and a perfect 
nesting condition are present, (a), (c), (e) and (g) show the 
calculated density of states whilst (b), (d), (f) and (h) are 
the calculated graphs of magnetisation as a function of I, the 
Stoner exchange parameter for various levels of band-filling, 
below half-filling. The tight-binding parameters for our dis- 
persion, e(k) = — cosfc^Q + p cos k y ) — 7 cos 2^ — ScosSkx 
are indicated. As in Fig. H, two magnetic transitions are often 
visible. 



to saturation magnetisation. According to Eq. il a tran- 
sition to saturation magnetisation would naturally kill 
any magnetically-mediated superconductivity due to the 
lack of one spin species. Thus we turn our attention to 
an upper transition not of saturating nature, and which 
could be first order, or close to first order, as the latter 
would allow the Fermi surface sheet configurations to be 
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FIG. 5: The phase diagram of our model in N, I space shows 
a line of first order transitions at I X (N) which end in a critical 
point. The transitions at the Curie onset, I c are also first or- 
der. This can be seen in the behaviour of (a) the longitudinal 
susceptibility and (b) the resulting (approximate) phase dia- 
gram. The tight-binding parameters used here are f3 — 0.7, 
7 = 0.03 and 8 — —0.03, in line with what follows. 



explored as the magnetisation is tuned by temperature 
or pressure, allowing the presence of sharp peaks in the 
density of states to be observed in measurements such as 
heat capacity and resistivity. 

The above considerations make the second example in 
Fig. ^ a most useful candidate, as we have the option 
of the upper transition being first order or close to first 
order, depending on N, and the magnetisation does not 
saturate at the upper transition. Henceforth, we take 
(3 = 0.7, 7 = 0.03 and 5 = -0.03. 

Interestingly, the magnetic phase diagram in the N, I 
plane seems to contain a line of first order transitions 
corresponding to I X (N). This line terminates at a crit- 
ical point at a value of N of about 0.76. For values of 
N below this, the transition in M(7) is smooth and can 
be described as close to first order. This is perhaps best 
seen in Figure |^(a), where we show the longitudinal sus- 
ceptibility, x as a function of I for (3 ~ 0.7, 7 = 0.03 
and 8 = —0.03, taking different values of N. We take the 
longitudinal susceptibility from the Stoner forms: 



Xfe 



p(pef) 
1 - Ip(/ief) 

2pOt)p(ah) 



Xpara 



p(p-[) + p(pi) ~ 2ip(pi)p(pi) 



(10) 

(11) 



in the paramagnetic (Fermi level hef) and ferromag- 
netic states, respectively. As seen in Figure ||(a), the 
susceptibility at I x grows as we move closer to the criti- 
cal point; the resulting (approximate) phase diagram in 
(N, I) space is shown in Figure |^(b). 



B. Magnetic field effects 

It has also been found that the features associated with 
T x and T c can be recovered at pressures above p x and p c 
respectively by the application of a magnetic field. In 
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FIG. 6: Our model predicts that both the magnetic jump at 
the onset of ferromagnetism and that at finite magnetisation 
can be pushed to lower values of exchange energy I c or I x 
respectively (corresponding to pressures higher than p c and 
p x ) by the application of a magnetic field, H. The tight- 
binding parameters used here are again [3 = 0.7, 7 = 0.03 
and S — —0.03. The number of spins, N is 0.8 in both figures 
(a) and (c) which show the calculated M(I) curves and the 
resulting phase diagram in H, I space. Recent experimental 
data, from Pfleiderer and Huxleyu is shown in (b) and (d). 
This has a form similar to that predicted in (a) and (c), with 
the appearance of two transitions in M, which are tunable 
with pressure. 



our model, this metamagnetism arises as a direct con- 
sequence of having spin-split Fermi surfaces, with the 
spin species being governed by number conservation and 
the free energy of the form given in Eq. [j]. As shown in 
Fig. ||, turning on the magnetic field, H, pushes both the 
magnetisation jump at I c (the Curie transition) and at I x 
(within the ferromagnetic state) to lower values of I, or 
equivalently, higher pressures. Both the predicted M(I) 
curves and the resulting phase diagram in H, I space, bear 
a striking resemblance to recent experimental datatl. We 
should point out that the calculated H, I phase diagram 
was obtained by looking for the maximum in the gradi- 
ent of the magnetisation for each M(I) plot in field. The 
first order nature of both transitions is softened with de- 
creasing I — corresponding to going to higher pressures. 
Experimentally, such softening is not observed at pres- 
sures reached this far, as can be seen from Figure g(b). 
We indicate where the magnetic transition in our model 
is no longer first order by a dotted line on the H, I phase 
diagram. 
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C. Superconductivity and the zero-temperature 
phase diagram 

From the heavy nature of the quasiparticles alone, we 
expect strong coupling theory to be required to gain a 
true estimate of Tsc in UGe2- However, we can gain use- 
ful information on the form of the phase diagram from 
examining the zero temperature, weak-coupling proper- 
ties of our model. We note here that, without a reli- 
able model of M(T), it is impossible to construct the 
usual Tsc equation in the BCS theory in the ferromag- 
netic state, as here we have an example of a temperature- 
dependent potential — the changing magnetisation alters 
the spin-dependent interaction. 

A reasonable estimate of the strength of triplet pairing 
can be gained from examining the relative magnitudes of 
the mass renormalisation parameter Xz and interaction 
parameter Aa, defined as 



Az = 



J FS J FS , <fk<Pk' Vr T (k-k') 

Ifs d 2 k 

_ J FS J FS , ^fcriWntk-kXkMk') 



, (12) 
, (13) 



where each integration is over the Fermi surface (FS) ei- 
ther in k or k' space. The term rj(k) is the angular part 
of the superconducting order parameter. In the above 
we are following the notation adopted by Monthoux and 
LonzarichEj and restrict ourselves to examine majority 
spin triplet pairing, using Vj-[-(q). The choice of order 
parameter should naturally reflect the symmetry proper- 
ties of the UGe2 crystal structure. Such considerations 
should lead us to examine non-unitary statesEZrEa, but 
here for simplicity we consider as an example the states 
Ak = Aq sm(k x ) and Ao sin(fcy), anticipating that which 
state is favoured may change as the majority Fermi sur- 
face is tuned through p x , making the transition from be- 
ing open to being closed. The favoured state is deter- 
mined by which i] (k) gives a larger positive value of Aa- 

There is a caviat at this stage. The sharpness of the 
van Hove singularities is useful in provid ing the magneti- 
sation behaviour shown in Section III A. However, as can 
be readily seen from Eq. |5|, it can also lead to a switching 
of the sign of the pair interaction from attraction to re- 
pulsion at zero q if the density of states at the van Hove 
point is too high. This is because Xcr^Cq = 0) = p(/j, a ). 
With the presence of a strict low dimensional van Hove 
singularity (an infinity in the DOS) at zero temperature, 
the quantity p{^)p(p\,) will be greater than 1 when the 
majority Fermi level reaches the van Hove point. 

There are several possible routes to softening the van 
Hove singularity. One would be to introduce a small 
amount of disorder into the model. Another would be 
to transfer the calculation to three electronic dimen- 
sions, rather than the current two, but that is com- 
putationally time-consuming. Here we take a simpler 



approach, calculating all XoS(q) at a small finite tem- 
perature (ksT = 0.045), retaining the original zero- 
temperature value of I for each magnetisation examined. 
This procedure maintains the order of the magnetic tran- 
sition at I x , although it does mean that the effective den- 
sity of states, Xo°t7 (q = 0) is inconsistent with that used 
to find M(I). Xctct (q) is calculated on a 40 by 40 grid for 
k and k' points in the Brillouin zone. The agreement is 
close away from the van Hove regions. 

In order to have two first order magnetic phase transi- 
tions (see Fig. ||(b)), we take N — 0.77 in what follows, 
and examine majority sheet, spin triplet superconductiv- 
ity. The final refinement to the model comes in the form 
of adjusting the q-dependence of I at this stage. Until 
now, we have considered a q-independent Stoner factor 
which arises from on-site repulsion of like-spin electrons. 
However, / should really convey some of the physics of 
electron-electron interactions at finite distances, and so, 
in calculations of the superconducting potential, we con- 
vert to the following form: 



(14) 



which effectively reduces high-q modes in the system, in 
line with the ferromagnetism of the real compound. This 
is the simplest, first approximation of such effects. 

In Fig. |?] we show a measure of the superconducting 
interaction strength, , for various values of £, our 

Stoner 'structure factor'. The term 1 + Xz measures the 
mass renormalisation — this is shown separately in Fig. 
||. Of course, jxf— is not the full story when one con- 
siders how the superconducting transition temperature 
Tsc should behave as a function of /. Roughly speaking, 

Tsc ~ w c e a a (see, for example, Fay and AppclO) 
and there is substantial variation of u> c , the paramagnon 
energy with /, especially when soft magnetic modes are 
present. The paramagnon energy can be expanded as 



cq 



(15) 



where C and c are temperature dependent parameters, 
equivalent to 7 and c in the work of Lonzaricbo. As de- 
veloped by Brinkman and EngelsbergcJ and implemented 
by Fay and Appel, the oj c prefactor will substantially re- 
duce superconductivity at the ferromagnetic I c and the 
same argument should apply around the secondary tran- 
sition at I x , although in both cases this effect is more 
drastic for a second order magnetic transition (where the 
susceptibility diverges, i.e. x ~^ 0) than for the first 
order cases here. What remains important, therefore, 
is the stable region of superconductivity in the ferromag- 
netic state, where jt^— is approximately flat and high, 
especially for higher Stoner structure factor, £. 

We might reasonably ask what it is that enhances su- 
perconductivity in the region between I x and I c . In this 
range of /, the majority Fermi level sits between the two 
van Hove peaks, and an examination of the q dependence 
of the interaction potential, Eq. || reveals that there is a 
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paramagnetic 




FIG. 7: A measure of the strength of superconductivity: 
as a function of Stoner interaction strength, 7, normalised 
with respect to 7 C , the value of 7 at the zero temperature Curie 
point. I x , the value of I for the second jump in magnetisation, 
akin to the pressure identified as p x in UGe2, gives rise to 
the peak at 7a:/7c ~ 1.34. Also shown is the magnetisation, 
scaled down by a factor of two in the dimensionless units we 
are using. The inset shows the ratio ^ is always less than 1, 
as expected. In both plots, we show the results for different 
values of 'Stoner structure factor', £ (Eq. |l^). The tight- 
binding parameters used here are again /3 = 0.7, 7 = 0.03 
and 5 = —0.03, and the number of spins, N = 0.77. 



paramagnetic 




5=0.02 
£•=0.03 
£•=0.05 

5=0.1 



FIG. 8: The mass renormalisation: 1 + \z as a function of 
Stoner interaction strength, 7, normalised with respect to 7 C , 
the value of 7 at the zero temperature Curie point. This quan- 
tity peaks in the ferromagnetic state and remains high with 
decreasing 7 into the paramagnetic state. Again, we show 
the results for different values of 'Stoner structure factor', £. 
The tight-binding parameters used here are again f3 = 0.7, 
7 = 0.03 and S = -0.03. The number of spins, N = 0.77. 



region of high, sustained zero q pairing for I c < I < 
At both ic and I x , V(q) is broad and high around q = 0. 
Outwith the region I c < I < I x the pairing potential is 
smaller and more localised around zero q. Furthermore, 
in the region between I c and I x , the mass renormalisa- 
tion, represented by Xz is also approximately flat and 
high. This mass enhancement, shown in Fig. M com- 
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FIG. 9: Evolution of majority (solid line) and minority (dot- 
ted line) Fermi surfaces as a function of Stoner exchange. 
We see the transition in the majority sheet from an open 
to a closed morphology, going through a van Hove singu- 
larity at 7 = I x . Also shown is the order parameter used, 
which changes from sin(fc s ) for 7 < I x to sixi{k x ) for 7 > I x . 
The tight-binding parameters used here are again /3 = 0.7, 
7 = 0.03 and S = —0.03, and the number of spins, N = 0.77. 



pares well with the high effective mass plateau found in 
de Haas van Alphen measurements on the ferromagnetic 
state between pressures p c and p x . 

Lastly we note that the order parameter, r)(k) favoured 
on tuning /, changes as we go through I x . This makes 
sense from the point of view of the Fermi surface topol- 
ogy, as depicted in Figure ||. 7y(k) = sin(fc y ) is favoured 
when the majority spin sheet is open — for all I < I x — 
and the sm(k x ) order parameter becomes favoured for the 
closed surfaces for I > I x . As already stated, these are 
overly simplified gap functions, so this change of symme- 
try may well be more gradual or masked in a non-unitary 
representation. 



IV. CONCLUSIONS 

In this paper we have proposed that the unusual phase 
diagram of UGe2 is a result of a novel tuning of the Fermi 
surface topology by the magnetisation of the ferromag- 
netic state. We have constructed a model for this which 
illustrates how superconductivity, the tunable magneti- 
zation features and the quasiparticle mass are related by 
a twin-peak density of states feature, consistent with ex- 
periment. Indeed, the twin-peak DOS feature is also ap- 
parent in the most-decent LDA +U calculations of the 
non-magnetic statea. 

We have first looked for an explanation of the two tran- 
sitions in the low temperature magnetisation of UGe2. 
We have shown how a density of states with two peaks 
as a function of energy will give rise generically to a 
zero temperature magnetisation which has two transi- 
tions as the ratio of Stoner exchange energy to bandwidth 
is tuned. We have demonstrated how this interesting 
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form of DOS arises when the bandstructure is quasi-one- 
dimensional, as is believed to be the case for UGe 2 . 

We are thus able to propose a mechanism for the 
asymmetry of the superconducting T$c with respect to 
the magnetic quantum phase transition associated with 
the Curie temperature, shifting the focus of our atten- 
tion away from p c towards p x as the magnetic transi- 
tion point of interest with regard to superconductivity. 
In our model, a large density of states at the majority 
spin Fermi surface is associated with the T x transition 
and is also driving superconductivity. This causes super- 
conductivity to be favoured in the ferromagnetic state 
relative to the paramagnetic state in a natural and pre- 
viously overlooked way, and could also be a reason for 
the enhancement of superconductivity in the ferromag- 
netic state of other materials. We further note that, in a 
real system, the presence of impurities will probably en- 



hance further the relative tendency to formation of the 
superconducting state in the ferromagnetic phase. 

There is also a straightforward explanation in the band 
magnetism scenario for the observed metamagnetic be- 
haviour and a large, constant mass enhancement between 
I x and 7 C , consistent in form with de Haas van Alphen 
measurements. Potentially we may have a reason for the 
disappearance of dHvA signal between p x and p c , which 
would here be due to the Fermi surface becoming open. 
However, we should note that only one set of authors 
have observed such disappearance of dHvA signalEI. 
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P.B. Littlewood, A.J. Millis, P. Monthoux, C. Pfleiderer, 
S.S. Saxena, T. Schmidt, M. Sigrist and T. Terashima 
for useful discussions. Financial support was supplied by 
the EPSRC. 
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By coupling high wavevector fluctuations to the q = 
mode, Watanabe and Miyake's model is reminiscent of 
the spin-bag mechanism of nodeless superconductivity pp 
posed by Schrieffer et al. in the context of the cuprates 
It is worth noting that, on scanning through energy in 
gMasi-one-dimensions, optimal nesting (which is also the 
density of states minimum) does not coincide with the van 
Hove point (the density of states maximum). This is be- 
cause optimal nesting in quasi-lD results from a Fermi sur- 
face consisting of parallel lines, each parallel to one axis 
of the Brillouin zone, as shown by the blue lines in Fig- 
ure ^(c). By contrast, the optimal nesting Fermi surface 
of a two-dimensional square lattice in a nearest-neighbour 
tight-binding model is a diamond shape, with parallel sur- 
faces at 45 degrees to the Brillouin zone axes. As shown in 
Figure |^( a), this Fermi surface also corresponds to the van 
Hove point and density of states maximum. 
Setting a x to 1 in Eq. ^ and varying I is equivalent to 
varying the quantity /' = I/a x , as can be seen from the 
form of the free energy in Eq. 
pressure is to increase, probab ; 
ping parameter a x . This suppresses ferromagnetism, as it 
reduces We show plots of M(I) since 1 = 1' here. 



j, as can De seen irom me 
ll. The effect of increasing 
y most strongly, the hop- 



